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ABSTRACT 

Infrared- Faint Radio Sources (IFRSs) are a class of radio objects found in the Australia Telescope 
Large Area Survey (ATLAS) which have no observable mid-infrared counterpart in the Spitzer Wide- 
area Infrared Extragalactic (SWIRE) survey. The extended Chandra Deep Field South now has even 
deeper Spitzer imaging (3.6 to 70 /im) from a number of Legacy surveys. We report the detections of 
two IFRS sources in IRAC images. The non-detection of two other IFRSs allows us to constrain the 
source type. Detailed modeling of the SED of these objects shows that they are consistent with high 
redshift (z > 1) AGN. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: starburst 



1. INTRODUCTION 

Infrared-Faint Radio Sources (IFRSs) were recently 
discovered in the Australia Te l escope Large Area Sur- 
vey (ATLAS) by iNorris et al] (|2006D . These are ra- 
dio sources brighter than a few hundred /iJy at 1.4 
GHz which have no observable infrared counterpart 
in the Spitzer Wide- area Infrared Extragalactic Sur- 
vey (SWIRE) (jLonsdale et al.|[200l . They may be re- 
lated to the o ptical l y invi sible radio sources found by 
iHigdon et al.l (|2005l . |2008( ). which are compact radio 
s ources with n o optical counterpart to i? ~ 25.7. 

INorri s et al] (.2006) and Middelberg et all (|2008a[) have 
identified 53 such sources out of 2002 radio sources in the 
ATLAS survey. Most have flux densities of a few hundred 
/xJy at 1.4 GHz, but some are as bright as 20 mJy. The 
IFRS sources were unexpected as SWIRE was thought 
to be deep enough to detect all radio sources in the local 
universe, regardless of whether star formation or Active 
Galactic Nuclei (AGN) powered the radio emission. Pos- 
sible explanations are that these sources are i) extremely 
high-redshift {z > 3) radio-loud AGN, ii) very obscured 
radio galaxies at more moderate redshifts (1 < z < 2), 
iii) lobes of nearby but unidentified radio galaxies, or iv) 
an unknown type of object. 

The nature of IFRSs has been hard to determine be- 
cause they have only been detected in the radio. Spec- 
troscopy is difficult because the hosts are optically faint 
and the radio positions can also have uncertainties on the 
order of a few arcseconds. Norris et al. (2006) stacked 
the positions of 22 IFRSs in the Spitzer 3.6 /im IRAC 
images and found no detection in the averaged image, 
showing that they are well below the SWIRE detection 



threshold. 

Rec ently, iMiddelberg et al.l (l2008bD and INorris et all 
()2007f ) targeted 6 IFRSs with the Australian Long Base- 
line Array (L BA) and successfull y detected two of the 
sources. The INorris et all (|2007f ) LBA detection con- 
strained the source size to less than 0.03 arcsec, sug- 
gesting a compa c t radio core powered by an AGN. 
IMiddelberg et al.l p)08b) found the size and radio lu- 
minosity of their LBA-dctccted source to be consistent 
with a high redshift (z > 1) Compact Steep Spectrum 
Source. The VLBI detections rule out the possibility 
that these particular IFRSs are simply the radio lobes 
of unidentified radio galaxies. iGarn fc Alexander! (|2008D 
stacked IFRS sources in the Spitzer First Look Survey 
at infrared wavelengths as well as at 610 MHz. They 
find that the IFRS sources can be modelled as compact 
Fanaroff Riley type II (FRII) radio galaxies at high red- 
shift (z ~ 4). Thus the evidence suggests that IFRSs are 
predominately high redshift radio-loud AGN. 

Ultra-deep Spitzer imaging is now available over the 30 
X 30 arcmin^ region of the Chandra Deep Field South 
SWIRE field. Four IFRSs lie in this region and in this pa- 
per we report on the constraints on the nature of IFRSs 
derived from the new Spitzer data. We assume a Hub- 
ble constant of 71 kms~^Mpc~^, and matter and cosmo- 
logical constant density parameters of VIm = 0.27 and 
^Ia = 0.73 in this paper. 

2. OBSERVATIONS, DATA AND SAMPLE SELECTION 

The extended Chandra Deep Field South is centered 
at 3h32m28s, -27° 48' 30". It overlaps the Great Ob- 
servatories Origins Deep Survey (GOODS) South field, 
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which is one of the best studied re gions of the sky. The 
ATLAS survey (|Norris et al.l 120061) consists of deep ra- 
dio observations of a 3.7 degree^ field surroundi ng the 
eCDF S which is coincident with SWIRE ( Lonsdale et aLl 
12004 . The ATLAS 1.4 GHz observations reach 20 ^iJy 
to 60 /iJy rms, with the deepest region covering the 30 x 
30 arcmin^ eCDFS. ATCA 2.4 GHz observations of the 
SWIRE region were obtained by the ATLAS team over 
the last two years and the current image reaches ~ 0.1 
mjy rms (Middelberg et al. in preparation). 

The SWIRE survey reaches 5cr depths of 3.5, 7.0, 
41, and 49 mJy , resp ectively, in the four IRAC bands 
(jLonsdale et al.ll2004i ). In the MIPS bands at 24 and 
70 fim the SWIRE 5cr depths are 189 iiJj and 16 mJy. 
To these depths, 22 radio sources in the full SWIRE 
Chandra Deep F ield South region were undetected, and 
classed as IFRS porris et al.ll2006f) . This paper focuses 
on four IFRS sources which lie in the 30 x 30 arcmin^ 
region of the extended Chandra Deep Field South. 

The eCDFS sub-region of SWIRE has been the target 
of two separate Spitzer Legacy proposals since SWIRE, 
and hence there are now deeper infrared data. The 
''Spitzer IRAC MUSYC Pubhc Legacy in E-CDFS (SIM- 
PLE) Legacy project (Damen et al. in prep) provides 
ultra-deep imaging of eCDFS in all four IRAC bands, 
and the "Far Infrared Deep Extragalactic Legacy (FI- 
DEL, PI: Dickinson) Legacy project adds ultra-deep 24 
and 70 /im observations. The 5cr depths of SIMPLE are 
about 0.8, 1.2, 6.3 and 6.6 fiJy in IRAG channels 1 to 4, 
respectively. FIDEL achieves 5a depths of 50 /iJy and 3 
mJy at 24 and 70/zm, respectively, in the areas with most 
coverage. We searched for counterparts to the four IFRS 
sources in the new Spitzer data. A matching radius of 
up to 3 arcsec was used for the IRAC and MIPS 24 fim 
data, and 8 arcsec for MIPS 70 fj,m. All four positions 
were examined by eye in all Spitzer bands for any obvi- 
ous failures (see Figure 1). We find two of the sources 
(S446 and S506) were detected in the SIMPLE ultra-deep 
IRAG imaging (see Figure 1) and we have improved IR 
constraints for the rest. Table 1 summarises the data for 
the four IFRSs. 

We have also searched for counterparts in existing 
deep optical a nd near-infrared (NIR) data of t he eCDFS 
from GOODS (Giavalisco et al.''2004'), GEMS ()Rix et al.l 
|2004) and MUSYG (Gawiscr ct al. 2006), and find an op- 
tical/NIR detection of only one IFRS source (S506). The 
GOODS AGS imaging reaches approximately Bab = 
28.1 , Vab = 28.9 , Iab = 28.3 and zab = 27 A. GEMS 
comprises two band ACS imaging and reaches depths of 
Vab = 28.25 and zab = 27.1. The MUSYC survey has a 
5(7 sensitivity of ab out Jab ^ 22.5 and Kab ^ 22.0 
(|Tavlor et al.i[2?)09l ). Two IFRSs, S415 and S506, lie 
within the deep optical/NIR imaging area and only one, 
S506, is detected as a very faint source. The lack of a 
detection in these optical and near-IR images for S415 
strongly implies that it is at high redshift and / or is very 
obscured. 

A 2 Ms X-Ray exposure of GOODS Sou th was obtained 
with the Chandra X -ray Observatory (iGiacconi et al.l 
I2002t iLuo et al.l 120081 ). The on- axis full band sensitiv- 
ity reaches 7.1 x 10"^'' erg s~^ cm~^ at the average aim 
point and the minimum full band sensitivity is about 3.3 
X lO^^^erg s~^ cm~^ over the GOODS-S region. Only 



one of the IFRS sources (S415) lies within the X-Ray 
coverage and there is no detection. 

The photometry of the IFRSs was fitted with the fol- 
lowing te mplates: 1) a p rototypical starburst ULIRG 
(Arp220, ' Silva et al.l[T99l . 2) a hot dusty AGN dom- 
inated ULIRG (Mrk231, phot ometry from NEEB), 3) a 
starforming M82-like galaxy (jCharv fc Elbad 120011 ) . 4) 
a radio-loud galaxy with an extreme radio-optical ratio 
(3C273, photom etry from NED) , and 5) an old stellar 
population from iMarastonI ()2005[ ). The old stellar tem- 
plates were fitted only to the two IF RS sources with 
IRAG detections (S446 and S506). The lMarasto^ (l20Q5l ) 
templates explored have a metallicity of 0.67, assume a 
Salpeter IMF, and have ages of 0.3, 1 and 5 Gyr. The re- 
sults of the Spectral Energy Distribution (SED) analysis 
are described in the next section. 

3. DISCUSSION 
3.1. S283 

This 0.287 mJy radio source lies just outside of the 
SIMPLE ultra-deep IRAG imaging and FIDEL 70 fiui 
imaging. There are no SWIRE optical data either, but 
the shallower SWIRE IRAC 3.6 and 4.5 /im imaging 
shows the radio source is positioned between two IRAC 
sources (Figure 1), which are 3.5 and 3.7 arcsec offset, 
respectively. The offsets are large (>5 times the uncer- 
tainty of the radio position) so we do not consider either 
IRAC source to be a plausible counterpart. The FIDEL 
24 /im imaging shows some possible fiux at the source 
position, but it is likely that this faint flux is confusion 
from the two nearby IRAC sources. To verify that there 
is no detection of the radio source in the IRAC and MIPS 
24 fj,m bands the two point sources were subtracted from 
the images and the residual image is shown in Figure [2j 
The residual flux at the radio position is 2.4 ± 1.8 /iJy 
at 3.6 /im and 13 ± 7 /iJy at 24 /im. This is consistent 
with the conservative upper limits we have assumed of 
3.5 /xJy at 3.6 /xm and 100 /iJy at 24 /xm. 

If this source was a typical ULIRG like Arp220 or Mrk 
231, even for ^ > 2, it would be detected by IRAC imag- 
ing (see Figure 3). For z < 2 the 24 /im flux density 
expected is >> 100 /iJy so it would be clearly detected 
in MIPS. A starforming galaxy such as M82 would also 
have clearly been detected by Spitzer. While a redshift 
less than 1 is allowed by the 3C273 SED, that would 
imply a low mass object. The stellar mass estimate at 
z = 1 is M ~ 6 X 10^ A/0 assuming the K band lumi- 
nosity fro m the 3C273 SED fit and a mass to light ratio 
of 1 (e.g. iBell fc de Jong! [200 It ). This is comparable to 
the most massive dwarf galaxies, which are not known 
to be strong radio emitters or contain AGN. The IRAC 
limits alone suggest a stellar mass of M < 6 x lO^^M© 
at z = 1. Radio-loud AGN with SEDs similar to 3C273 
are found in massive galaxies with M > 10" M©. The 
IRAG non-detection therefore suggests S283 is a radio- 
loud AGN at z > 1. This is supported by the galaxy's 
MIR-radio correlation, q = log(S'24^m/S'i.4GHz) < —0.46, 
whic h is well in the ra nge expected for radio-loud AGN 
(e.g. IBovle et al.l[20?)7l ). 

1 http://nedwww.ipac.caltech.edu/. The NASA/IPAC Extra- 
galactic Database (NED) is operated by the Jet Propulsion Labo- 
ratory, California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 
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3.2. S415 

This IFRS has a flux density of 1.21 mJy at 1.4 GHz, 
but is undetected in all Spitzer bands. Since it lies in 
the GOODS-S proper, where the data is most sensi- 
tive, it has by far the most extreme flux density ratios 
of the sources discussed here. It is as extreme, if not 
more s o, than the optically faint but submm bright HDF 
850.1 (jPunlop et al.l[200l iCowie et all 120091 ). This ra- 
dio source has no ACS counterpart and we measured the 
3cr ACS limits at the source position to be Bab ^ 28.1. 
Vab > 28.9, Iab > 28.3, and zab > 27.4, using a 0.6" 
diameter aperture. There is no near-infrared counterpart 
in the MUSYC or VLT ISAAC observations, giving 3cr 
limits of Jab > 25.5, Hab > 25.8 and Ksab > 25.5. 

If this source has an SED similar to M82, Arp220 or 
Mrk 231 it would be detected in all Spitzer bands, even 70 
jum, for z < 4. If the radio-loud 3C273 SED is adopted, 
this source must lie at z > 3 for it to remain undetected 
in the IRAC and ACS bands (assuming no obscuration, 
see Figure 4). At redshift z = 1, the IRAC 3.6 and 4.5 
/im limits are a factor of 2 to 3 below the 3C273 SED, 
which implies an extinction of 0.8 to 1.1 mags would 
be required for this source to be undetected in the IRAC 
channels. This level of extinction i s seen in some extrem e 
ULIRGs (e.g. IGenzel et all [T998HM^hv et al.l[200l . 
but this source is not luminous in the infrared. This 
source therefore lies at redshift z >> 1. Even at these 
redshifts obscuration of a few magnitudes is required for 
a 3C273 object to be undetectable in the optical/NIR 
bands (see Figure 4). At z > 1 the radio luminosity of 
S415 is Pi.4GHz > 5 X 10^^ W Hz-^ The Spitzer non- 
detection of S415 therefore suggests that this source is a 
distant (z >> 1) obscured radio- loud AGN. 

The X-Ray non-detection implies that this source has 
an X-Ray luminosity Lo.5-8fcey ^ 2 x 10^^ erg at 
z = 2 and Lp s -skeV ^ 6 x 10^^ erg s^^ at z = 3, where 
the lLuo et all (2008') limit assumes an X-Ray power law 
of r = 1.4. This suggests that, for these redshifts, either 
the source has a very high column density of neutral hy- 
drogen obscuring the X-Ray emission, which is consistent 
with the optical/NIR non-detection, or the source is in- 
trinsically X-Ray faint for an AGN. 

A preliminary analysis of ATCA 2.4 GHz data (Mid- 
delberg et al. in preparation) finds S415 has a significant 
detection (> 7a) of 0.67 mJy. This implies a steep ra- 
dio spectral slope of a — —1.10, which can not be pro- 
duced by star forming processes. S415 has a slope similar 
to some ultra steep spectru m sources (Rottgcrin g et al] 
Il997t iDe Breuck et all 120041 ) . which have been linked to 
massive high-redshift (z > 2) radio galaxies. 

3.3. S446 

IFRS S446 is a 0.338 mJy radio source which lies just 
outside the GEMS ACS coverage and MUSYC near- 
infrared imaging. There is a faint source in the IRAC 
3.6 and 4.5 /im channels just 2.2 arcsec north of the ra- 
dio position (see Figure 1) and we assume this is the 
counterpart to S446. The probability that one or more 
IRAC sources lies randomly within a distance 6* of a ra- 
dio source is P = 1 — exp(— 7m0^), for an IRAC source 
density n (often called the P-statistic; e.g. iDownes et al.l 
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I1986D . The P-statistic is only a rough estimate because 
it does not does not take into account the individual 
positional uncertainties and assumes a random distribu- 
tion of the background population, whereas astronomi- 
cal sources are clustered. Nevertheless, for the SIMPLE 
IRAC source density of 116700/deg^ the P-statistic sug- 
gests that the chances of an IRAC source lying closer 
than 2.2 arcsec is 13%, so the counterpart is reasonably 
reliable. The source is 6.6 ± 0.3 /iJy and 5.7 ± 0.5 ^Jy 
at 3.6 and 4.5 /im, respectively, and not detected in the 
other IRAC bands. There is no detection in the 70 and 
24 fim images, but limits are hard to quantify because 
the source falls in between 2 brighter sources. We assume 
the noise here is twice that of the local noise in the MIPS 
imaging for the purposes of the SED fitting. 

Similar to the other IFRS sources, the non-detection 
in the two longer wavelength IRAC bands, and MIPS 
24 and 70 fiui bands, rules out M82, Arp220 and Mrk 
231 SEDs for this source, at a wide range of redshift 
(z < 6). A 3C273 SED would not produce the MIR 
emission detected in the IRAC bands (see Figure 5), so a 
possible explanation for this source is a radio-loud AGN 
which dominates the radio emission but with a stellar 
compo nent which is see n in the MIR. 

The iMarastonI (|2005[ ) old stellar population fit to the 
IRAC data constrains the redshift of this galaxy to 1 
- 1.5, and the best fit is a 1 Gyr old model at z = 1.5. 
At this redshift the radio luminosity of S446 is Pi.4ghz = 
3.7 X lO^'' W Hz^^, which would place it at the low power 
end of r adio-loud AGN. 

Using IMarastonI ()2005[) SEDs with a different metal- 
licity, including a diff erent library of SEDs (e.g. 
iBruzual fc Chariot! |2003() . or adding reddening as a free 
parameter would give a different best- fit redshift, but this 
level of complexity can not be explored with only two 
datapoints. Lastly we note that the IRAC detections are 
unlikely to be caused by the hot dust component of an 
AGN torus. The IRAC detections imply the MIR peak is 
shorter than 3.6 /im, so the dust temperatures would be 
greater than ^1600 K at z = 1 and greater than ~2400 
K at z = 2, using Wien's law. Silicate grains sublimate 
at about 1000 K, while graphite grains sublimate around 
1500 K, so these temperatures are too high for an AGN 
torus. Instead, AGN torus models typi cally show a cooler 
MIR emission peak of 7 to 10 /im (e.g. iSchartmann et aLl 
120051) . which is longward of IRAC channels 1 and 2. 

3.4. S506 

This source is similar to S446. This 0.170 mJy ra- 
dio source has a GEMS ACS V and z band detection 
1.4" south of the radio position that is 26.27 and 25.62 
AB magnitudes, respectively (see Figure 1). The likely 
counterpart in the IRAC 3.6 and 4.5 /im channels is less 
than 1" away from the ACS source, roughly 2.3" south 
of the radio position. This ACS counterpart, while faint, 
is a 7.0 and 5.1 sigma detection in the V and z bands, 
respectively. 

The P-statistic suggests the chance of this IRAC 
source being a random coincidence is about 14%. The 
source is 5.5 ± 0.3 /iJy and 5.5 ± 0.4 fiJy at 3.6 and 4.5 
/im, respectively, and not detected in the other IRAC 
bands. There is no MIPS 70 or 24 /im detection in the 
FIDEL images of this source. The FIDEL image shows a 
possible 24 /xm excess but this has a flux density of only 
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7.3 ± 4.2 /iJy and it is not coincident with the IRAC 
source. There is also no detection in the MUSYC NIR 
imaging {Jab > 22.5 and Kab > 22.0). Again, M82, 
Arp220 and Mrk 231 type SEDs are ruled out for this 
source by the non-detection in the MIPS 24 and 70 fim 
bands. The MIR peak is between 3.6 and 4.5 fj,m and 
if hot dust is responsible for this then Wien's law im- 
plies hot dust temperatures of '^1400 K at redshift 1 and 
-2200 K at redshift 2. Similar to S446, this rules out the 
IRAC detection of a hot AGN tori. The most likely ex- 
planation for this source is a radio-loud AGN at (z > 2 ) 
with a stellar component dominating the IRAC channels 
(Figure 6). 

Using t h e IR AC data alone, we find the best fit 
iMarastonI (|2005i ) old stellar population model to the 
IRAC data is 1 Gyr old and places this galaxy at z — 2.5. 
The 5 Gyr model has a redder SED that would be de- 
tected by IRAC 5.8 /im imaging at z < 2, and the 1 Gyr 
model becomes a bad fit for z > 2.6. However the ACS 
detection suggests a blue excess from star formation ac- 
tivity. The best fit stellar popula tion model to bo th ACS 
and IRAC detections is a 0.3 GyrlMarastonl (pOOl model 
at z = 2.0. As for S446, the small number of datapoints 
does not warrant exploring all the free parameters in the 
redshift fitting, such as metallicity and reddening. How- 
ever, SEDs younger than 0.3 Gyr were also explored for 
S506 and these young SEDs have a blue excess inconsis- 
tent with the faintncss of S506 in the optical bands. At 
the best fit redshift of z = 2.0, S506 has a radio luminos- 
ity Pi.4GHz = 3.6 X 1024 W Hz-\ which, similar to S446, 
would place it at the low power end of radio-loud AGN. 

The ACS counterpart is made up of two clumps (see 
Figure 1), one of which is extended and has faint emission 
over 0.3 arcsec, which is 2.5 kpc at redshift z = 2. The 
dual nature of S506 could be explained by one component 
having an AGN while one (or both) has some recent star 
formation. 

3.5. IR-radio Correlation 



The IR-radio correlation (e.g. lYun et al.|[200TI ) IS an in- 
dicator of the dominant emission mechanism in a galaxy 
because both IR and radio emission are thought to be 
strongly linked to star formation. Any deviation from 
the correlation is a sign of an AGN. Galaxies with excess 
radio emission probably contain a radio-loud AGN, while 
IR-excess sources are likely to be radio-quiet AGN with 
hot dust dominating in the MIR (for z > 1). The ob- 
served MIPS to radio flux density ratio limits are shown 
in Figure 7. Some of the xFLS galaxies (Frayer et al. 
2006) have ratios consistent with the IFRSs but all are 
easily detected at optical and infrared wavelengths. The 
'S'24/im/'S'i.4GHz limits for all 4 IFRSs show that they have 
excess radio emission for z < 2. The limits from 70 fim 
are not as stringent as that from 24 /im, but they do 
show that the IFRSs have as much excess radio emis- 
sion as the AGN-dominated ULIRG Mrk231. The ex- 
treme source S415 has a ratio that is consistent with a 
radio-loud AGN at z > 4. The IFRS IR-radio ratios 
are consiste nt with a sample of hi gh redshift radio galax- 
ies (HzRGs. ISevmour et al.|[2007[) . which host luminous 
radio-loud AGN. 

3.6. Could IFRSs be galactic objects? 



We now consider whether IFRSs may be galactic ob- 
jects such as pulsars or radio stars. The majority of 
known pulsars are young objects with a low spin rate 
and these are distributed across the Galactic disc, as ex- 
pected for objects that have originated fairly recently 
from massive stellar supernovae. They are rare objects 
at high galactic latitude. For example, the Parkes Multi- 
beam Pulsar Survey, with a 1.4 GHz flux density limit 
of 0.15 mJy for long per iod pulsars, finds a density of 
about l/deg2 for |&| < 1 ()Camilo et al.ll2000l ). The den- 
sity of these pulsars drops to <0.25/deg2 by |5| = 4 deg. 
Applying the P-statistic from Section 3.3, we find the 
probability of one or more pulsars within 15 arcmin (half 
the eCDFS size) of a radio position is less than 5%, as- 
summg a pulsar densityof <0.25/deg2. The eCDFS has 
a galactic latitude of about -55 degrees, so the space den- 
sity of pulsars in the eCDFS field is much less than that. 
Therefore the 0.25 degree eCDFS field is unlikely to con- 
tain a pulsar. 

Shorter period pulsars, so-called millisecond pulsars, 
are thought to have been spun up to high rotation rates 
as a member of a low-mass X-ray binary (LMXB) system. 
These objects are usual ly found in globular clusters (e.g 
[Manchester et al.l 12001] ). The eCDFS is not a globular 
cluster field and so is unlikely to contain a LMXB. 

The optical non-detection of the radio sources gives 
us another clue about whether the sources are pulsars. 
Neutron stars have been detected in the optical, and 
we can take two examples: the Crab pulsar and the 
Vela pulsar. The Crab pulsar's optical counterpart is 
known as Baade's star and it is relatively bright in the 
optical, wi th V = 16.6. The Crab pulsar lies at 2 kpc 
fManchcst er et al] [2005'). so at 10 kpc it would have V = 
20.1 and V = 25.1 at 100 kpc. Thus, if these IFRSs are 
as bright as the Crab pulsar in the optical they would be 
easily detectable in the optical imag ing. The Vela pulsar 
is much fainter however, V = 23.6 (jMigna ni fc C aravec 



2001), and lies at a distance of 294 pc (Caraveo et al 
20011 ). An object such as this would be too faint for 



eCDFS HST detection if it lies further than about 2.3 
kpc. 

X-Ray emission from the neutron star surface or from 
the the pulsar wind nebula maybe detectable if the IFRS 
is a pulsar. The one source with X-Ray data, S415, has 
an X-Ray luminosity limit that is 100 times less than 
some of the faintest X-Ray emitting pulsars known (e.g. 
IKargaltsev fc Pavlov|[2009l ). X-Ray emission from a neu- 
tron star may not be orientated the same way as the 
radio emission, and the X-Ray-radio observations were 
not simultaneous, so this is not conclusive evidence that 
S415 is not a pulsar. It is however consistent with the 
idea that S415 is an extragalactic source. 

Could IFRSs be main sequence stars? Ultracool dwarf 
stars, a class thought to be radio active, were observed 
at 4.8 GHz, but only one was detected in a survey of 
eight ()Antonova et al.l l2008'). This star lies at a distance 
of 12.2 pc an d has a 4.8 GHz flux density of 0.286 mJy. 
iBergeiH ( 2006( ) observed 90 M and brown dwarf stars and 
found 8.5 GHz flux densities well below 1 mJy, although 
flares can increase the flux density by at least a factor of 
a few. The stars in this sample are all closer than about 
13 pc and have J and K magnitudes brighter than 18 
mag. The typical M-dwarf has an absolute magnitude of 
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8 - 17 My (e.g. lKaleiill997D . so at a distance of 1 kpc (10 
kpc) M-dwarfs have an observed magnitude of V — 18 
- 27 (23.1 - 32.1). Therefore the weakest of these would 
be undetected in the optical images at large distances, 
but the radio studies suggest M-dwarf stars close enough 
to be detected in the radio would be seen in the optical 
images. 

3.7. Comparison to other galaxy populations 

We can compare the space density of our sample of 
IFRSs with that of other high redshift z ~ 2 samples. For 
example, much work has been done on the BzK selection 
technique (jPaddi et al.ll2004 ) to select star- forming and 
passive galaxies at z ^2. The space density of BzKs 
is ~ 1 per arcmin (IDaddi et al.l 120071) . more than 200 
times greater than that of the IFRSs in this work (^16 
per deg^). 

Another high redshift sample is Dust-Obscured Galax- 
ies (DOGs), which are selected using a combination of red 
colors {R— [24] > 14, in Vega magnitud es) and bright fiii x 
densities in infrared (524 > 0.3 mJy) ()Dev et al.ll2008[ ). 
These criteria proved remarkably efficient for selecting z 
~ 2 galaxies. DOGs contribute ~ 26% of the total IR 
luminosity density at z = 2 and 60% of the total from 
ULIRGs. They have a space density of 0.1 per arcmin^ 
or 2.82 ± 0.05 x 10~^ h'^ Mpc~^, similar to submilhme- 
ter galaxies and about 20 times more numerous than the 
IFRSs. 

The paucity of IFRSs is not a surprise as they were 
selected to have extreme radio to infrared ratios, and as 
such are rare compared to other populations selected in 
the optical and infrared. How do AGN samples selected 
in the radio compare? By matching Faint Images of the 
Radio Sky at Twenty-cm (FIRST) 1.4 GHz r adio sources 
with the Sloan Digital Sky Survey (SDSS), llvezic et"al] 
(|2002D find about 3 radio selected AGN per deg^ with 
5*1.4 > 1 mJy and i < 21. The four IFRSs in this study, 
likely radio-loud AGN at z > 1, have a source density 
~16 per deg^. This suggests that IFRSs are a popula- 
tion of radio-loud AGN which have been unstudied by 
previous radio work. 

S446 and S506 lie at redshifts 1 to 2 and have radio 
powers (Pi.4GHz 

) of - 10^4 WHz^i, The loc al luminos- 
ity fu nction of radio AGN ([Best et al . 2005; S adler et al.l 
120021) suggests that there are 40 to 60 AGN per deg^ 
with radio powers of lO^'' - 10^^ W Hz^^ in the volume 
between redshifts 1 to 2, assuming n o evolution. If lu- 
mino sity evolution (e.g. IDonoso et al.|[2009 : Sadl er et al.l 
|2007[) is applied the number of AGN with radio pow- 
ers of 10^4 -10^5 W Hz-i increases to 110 to 140 per 
deg^. Therefore, the IFRSs presented in this work make 
up only a small portion (< 1%) of the total AGN pop- 
ulation with similar radio powers at that epoch. The 
remainder are presumably bright enough to be detected 
in the Spitzer bands. 

The IFRSs can also be compared to high redshift radio 
galaxies(HzRGs). A representative sample o f 69 HzRGs 
at z > 1 were observed with Spitzer (Sev mour et al.l 
I2007D and these massive galaxies (10^^ - 10^^'^ Mq) have 
the mid-IR luminosities of LIRGs or ULIRGs. While 
HzRGs have much greater radio flux densities than 
IFRSs, they have extreme MIR-radio ratios consistent 
with the available data on IFRSs (Figure 7). Hence, 



IFRSs maybe the lower luminosity analogs of HzRGs. 
There are only a few hundred known HzRGs across the 
full sky, so they are even rarer than IFRSs. 

4. SUMMARY 

We have searched for infrared counterparts to four 
Infrared Faint Radio Sources in the extended Chandra 
Deep Field South field using recently available ultra-deep 
Spitzer observations. No IFRS is detected in ultra-deep 
24 fim imaging, implying that the sources do not fol- 
low the IR-radio correlation and star formation can not 
produce all the radio emission observed in these objects. 
We find IRAC detections for two of the sources, and the 
non-detections of the other two provides constraints on 
the source SEDs. Typical ULIRG SEDs, such as Mrk231 
and Arp220, and galaxies are ruled out by the Spitzer 
data. The most likely explanation for these sources is 
that they are radio-loud AGN, with radio-to-optical ra- 
tios similar to 3C273, at redshifts z > 1. The most ex- 
treme source (S415) lies at z » 1 and requires several 
magnitudes of obscuration in the optical/NIR to remain 
undetected by deep imaging. It is very unlikely that 
the IFRSs are Galactic sources, but current data cannot 
conclusively rule this out in the case of the undetected 
sources. 

For the two sources with IRAC detections we find the 
SED can be described by a radio-loud 3C273-like SED 
combined with a stellar population. The stellar popu- 
lation SED fits to the IRAC MIR data, combined with 
the available optical limits, suggest that the two detected 
IFRS sources have redshifts of z ~ 1.5 and z ^ 2.0. For 
one source, S506, the ACS detection suggests the galaxy 
may have a stellar component that is only 0.3 Gyr in 
age. These two sources are similar to the IRAC detected 
OIRSs which are posited to b e "red and dead" radio 
galaxies at z > 1 ([Higdon et al.| [2b08V 

The evidence is mounting that a significant proportion, 
if not all, of the IFRS sources are radio-loud AGN at high 
redshift (z > 1), and not merely lobes of an unidentified 
radio galaxy. The sourc e density in eCDFS is ~1 6 per 
deg2 for 5i.4 > 0.1 mJy. IGarn fc Alexandeil (pOOSh find 
a source density of ^^3.5 per deg^ for 5i.4 > 0.5 mJy. 
While these sources are rare, they point to a population 
of AGN at high redshift that has been undiscovered until 
recently. 

This work is based in part on observations made with 
the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, California Institute of Tech- 
nology under a contract with NASA. Support for this 
work was provided by NASA through an award issued 
by JPL/Caltech. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) which is op- 
erated by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National 
Aeronautics and Space Administration. 
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S283 


S415 


S446 


S506 


RA (J2000) 


3;30;48.686 


3:32:13.077 


3:32:31.540 


3:33:11.486 


Dec (J2000) 


-27:44:45.32 


-27:43:51.07 


-28:04:33.53 


-28:03:19.09 


radio 1.4 GHz 


0.287 mJy 


1.21 mJy 


0.338 mJy 


0.170 mJy 


radio 2.4 GHz 


< 0.40 mJy 


0.67 mJy 


< 0.45 mJy 


< 0.45 mJy 




< 0.6 


-1.1 ± 0.13 


< 0.5 


< 1.8 


ACS B mag 




>28.1 






ACS V mag 




>28.9 




26.27 


ACS I mag 




>28.3 






ACS z mag 




>27.4 




25.62 


SWIRE g' mag 






>25.3 




SWIRE r' mag 






>24.8 




SWIRE i' mag 






>23.7 




J mag 




>25.5 




>22.5 


H mag 




>25.8 






K mag 




>25.5 




>22.0 


IRAC 3.6 /im 


< 3.5 /iJy 


< 0.8 /iJy 


6.6 ± 0.3 AtJy 


5.5 ± 0.3 /iJy 


IRAC 4.5 /im 


< 7.0 AiJy 


< 1.2 /iJy 


5.7 ± 0.5 AtJy 


5.5 ± 0.4 //Jy 


IRAC 5.8 /im 


< 41 /iJy 


< 6.3 /nJy 


< 6.3 (UJy 


< 6.3 /iJy 


IRAC 8.0 /im 


< 49 nJy 


< 6.6 /uJy 


< 6.6 /uJy 


< 6.6 /iJy 


MIPS 24 /im 


< 100 AtJy 


< 50 /uJy 


< 80 nJy 


< 80 //Jy 


MIPS 70 /im 


< 16 mJy 


< 3 mJy 


< 5 mJy 


< 5 mJy 



TABLE 1 

Summary of the observed properties of the four Infrared Faint Radio Sources in the extended Chandra Deep Field South. 



3.6 fim 



4.5 fim 



24 fim 



ACS V 



S283 



S415 



S446 



S506 




X 



. V ■ _ ■ ■ » / 



Fig. 1.— Grey scale Spitzcr and ACS images of S283, S415, S446 and S506. From left to right, images are IRAC 3.6, 4.5 ^im, MIPS 24 
/im, and ACS V band. The ACS grey scale is inverted for clarity. The green contours are the radio 1.4 GHz image with levels set at 50, 
100, 150, 200 fj,3j. The red cross marks the radio position. The Spitzer images are 25 X 25 arcsec and the ACS image is 12.5 X 12.5 arcsec. 
The IRAC counterpart to S446 can be seen 2.2 arcsec north of the radio position. The red circle in the S506 ACS image mark the very 
faint optical counterpart, which is consistent with the proposed IRAC counterpart 2 arcsec south of the radio position. 




Fig. 2. — The residual image for S283 with the two point sources near the radio center removed, for IRAC 3.6 ^im (left) and MIPS 24 
fim (right). As for Figure 1, The green contours are the radio 1.4 GHz image with levels set at 50, 100, 150, 200 fijj and the red cross 
marks the radio position. There is very little IRAC or MIPS flux at the radio position. 
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Fig. 3.— The SED of S283, fit by Arp220, Mrk231, M82 and 3C273. Points in black assume the galaxy is at z = 2, and blue are the 
data for 2: = 1. This galaxy is not well fit by M82, Arp220 or Mrk231, as it would be detected by IRAC and 24/^m imaging. The SED of 
radio-loud galaxy 3C273 at z > 1 is consistent with the observations. The SEDs are scaled up in luminosity by 90, 12 and 2000 times for 
Arp220, Mrk231 and M82 respectively, and 3C273 is scaled down by a factor of 1800. 
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Fig. 4.— The SED of S415, fit by Arp220, Mrk231, M82 and 3C273. Points in black assume the galaxy is at z = 4, and blue are the 
data for z = 3. This galaxy is not well fit by M82, Arp220 or Mrk231, as it would be detected by Spitzer imaging. The SED of radio-loud 
galaxy 3C273 at z >> 1 is consistent with the observations, if obscuration at optical wavelengths is assumed. The SEDs are scaled up in 
luminosity by 2.7 X 10'^, 3.9 X 10^ and 6.0 X lO'^ times for Arp220, Mrk231 and M82 respectively, and 3C273 is scaled down by a factor 
of 30. 
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rest wavelength (mlcronB) 
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Fig. 5.— Left: The SED of S446, fit by Arp220 Mrk231, M82 and 3C273. The SEDs are scaled up in luminosity by 750, 110 and 17000 
times for Arp220, Mrk231 and M82 respectively, and 3C273 is scaled down by a factor of 550. An old stellar population (red line) from 
[Maraston (2005) is fit to the IRAC data. Points in black assume the galaxy is at 2 = 4, and blue are the data for z = 1.5. This galaxy 
is not well fit by M82, Arp220 or Mrk231, as it would be detected by 24 fira imaging. The SED of radio loud galaxy 3C273 at z ~2.0 
combi ned with the old stellar population is consistent with the observations. Right: A zoom of the optical/NIR region. The IMarastonI 
II2005I ) stellar population models are shown at the best fit redshift of 1.5. The 0.3 and 5 Gyr models are satisfactory fits, but the 1 Gyr 
model provides the best fit. The optical limits arc g, r and i limits from SWIRE. 



Arp??0, z = -1 
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Fig. 6.— Left: The SED of S506, fit by Arp220, Mrk231, M82 and 3C273. The SEDs are scaled up in luminosity by 380, 55 and 8500 
times for Arp220 , Mrk231 and M82 respectively, and 3C273 is scaled down by a factor of 530. An old stellar population (red line) from 
IMarastonI I I2005I) is fit to the IRAG data. Points in black assume the galaxy is at z = 4, and blue are the data for z = 2.0. This galaxy 
is not well fit by M82, Arp220 or Mrk231, as it would be detected by 24 /xm imaging. The SED of radio loud galaxy 3C273 at z ^ 2.0 
combined with the old stellar population is consistent with the observations. Right: A zoom of the optical/NIR region. The IMarastonI 
l|2005) stellar population models are shown at the best fit redshift of 2.0. The 1 and 5 Gyr models fit the IRAG data, but only the 0.3 Gyr 
model (black line) can reproduce both IRAG and AGS detections. The J and K limits are from MUSYG. 
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Fig. 7. — The ratio of Spitzer 70 (left) and 24 (right) to 1.4 GHz flux densi ty, q, as a function of redshift. The four IFRSs are shown 
as black arrows at z = 0.06, as we have limits only. Red circles are xFLS galaxies lIFrayer et al."200(?), which ar e predominately star forming 
galaxies. Blue boxes or arrows show ratios from a sample of high redshift radio galaxies (Seymour et al.ll2007l ). which are r adio-loud AGN. 
The black SED tracks are Arp220 (dotted), Mrk231 (dashed) and M82 (dot dashed). Blue solid lines show tracks from lGharv fc Elbad 
1)20011 ) with bolometric IR luminosities ranging from 10** to 10^^ Lq . 



